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3.3A AN EVALUATION OF THE ACCURACY OF SOME RADAR WIND PROFILING TECHNIQUES 

A* J. Koscielny and R. J. Doviak 

National Severe Storms Laboratory 
1313 Halley Circle 
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INTRODUCTION 

Major advances in Doppler radar measurement in optically clear air have 
made it feasible to monitor radial velocities in the troposphere and lower 
stratosphere. For most applications we want to monitor the three dimensional 
wind vector rather than the radial velocity. Measurement of the wind vector 
with a single radar can be made assuming a spatially linear, time invariant wind 
field. The components and derivatives of the wind are estimated by the para- 
meters of a linear regression of the radial velocities on functions of their 
spatial locations. The accuracy of the wind measurement thus depends on the 
locations of the radial velocities. 

PETERSON and BALSLEY (1979) point out that a tradeoff exists for a given 
technique between the accuracies of horizontal and vertical component measure- 
ments. Because we usually need to measure the three components of wind with 
different accuracy and as inexpensively as possible, we are led to evaluate the 
suitability of some of the common retrieval techniques for simultaneous 
measurement of both the vertical and horizontal wind components. The techniques 
we will consider are fixed beam, azimuthal scanning (VAD) and elevation scanning 
( VED) . 

ERROR ANALYSIS THEORY 


The estimation of the parameters of a linear wind field from radial velo- 
cities is discussed by KOSCIELNY et al. (1982). The measured radial velocity v r 
can be modelled by a linear regression equation of the form 
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( 1 ) 


where P m is a row vector of regressor variables which are functions of range r, 
azimuth 9, and elevation angle 0 e » Km is a column vector of m parameters. The 
measured v r , a reflectivity weighted mean of radial velocities within the 
radar 1 s resolution volume, can contain errors £ due to nonuniform reflectivity, 
turbulence, targets such as hydrometeors that move relative to the wind, and a 
nonlinear wind. It can be shown that, given n measurements of v r , least 
squares estimates of K^ are computed by 
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where T indicates transpose and P^ is an nxm matrix of the regressor variables 
corresponding to the n radial velocity measurements in V n . Measurement errors 
in the radial velocities produce uncertainties in the estimate Km and the 
covariances of about K^ are given by 
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If the wind field has variations not modeled by (1), the will be biased 
and the amount of bias B m is given by the product of a known alias matrix 
with the vector of the % unknown parameters of the wind field not included 
in K m . Thus 
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where 


A mJl “ ^ P nm P nm^ ^ P nm P nP‘ 


(5) 


is a matrix of regressor variables for the components not included in (1). 

The various techniques referred to in the introduction assume the wind to 
be uniform (i.e. , the first and higher order derivatives are zero) over the data 
analysis volume. However, because w cannot be uniform for any appreciable depth 
of the troposphere (i.e., w must be zero at the earth’s surface), the horizontal 
wind can never be uniform at all heights. Thus we must account for errors 
produced by wind shear. We propose to analyze the errors in these techniques by 
computing the bias and variance of the least squares estimates with 
assumptions that is constant and the wind field is actually linear. Thus K~ 
contains the three uniform components u Q , v 0 , w Q and Kg contains the 8 spatial 
derivatives (u x , u z , v , v z , u y + V x , w x# w y , w z ) of the linear wind. In our 
evaluation and comparison of techniques, we assume that a total of n measure- 
ments are available for each and these n measurements are distributed in space 
to estimate wind at some height h. 

(a) Fixed Beam 

We consider a configuration for the fixed beam technique in which three 
beams, one verticl and two off-vertical at elevation 8 e , are sampled. The off- 
vertical beams usually have perpendicular horizontal projections; for 
convenience, we will consider them to have azimuths 0° and 90°. The total 
number of radial velocity measurements for a height h for all three beams is n; 
for generality, we let the number of vertical measurements be N. 

The bias and variance properties of the estimates K^- (u Q , v Q , w Q ) are 
computed in Appendix 1 using (3) and (5), and we find that, for n = 3N, 
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The bias equation is approximate because we have used h~rsin0 e which should be 
appropriate for r<30 km. From (7) we see that the bias due to spatial 
derivatives is a linear function of height which is expected because the beam 
separation is linearly dependent on h. In addition, we see from (6) and (7) 



194 


that the variance decreases with n, but that the bias cannot be reduced by data 
averaging. 

(b) Azimuthal Scanning 

In an azimuthal scanning technique, usually called VAD (Velocity 
Azimuth Display), data along a circle centered on the radar are used to 
directly estimate the components of the uniform wind field. The results of the 
bias and variance equation evaluation, in Appendix 2, are that 
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(c) Elevation Scanning 


In the elevation scanning (Velocity JElevation Display) technique, 
radial velocities at a height h are collected for elevation angles 
0 o £0 e £18O-8o* We assume data are collected for the two azimuths 0° and 90° 
so both horizontal components are measured. It is shown in Appendix 3 that 
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ERROR COMPARISON 


The results of our analysis of the three techniques, summarized in Table 1 
show the variances of the wind estimates all depend on a Z /n. Since the 
variance of an average of n independent data is a ^/n, w& will divide this 
quantity by the variance of the estimate of the wind component. Because of its 
similarity to the usual statistical definition, we term this quantity the 
efficiency of the estimate. 

The variation of the efficiences of the horizontal wind estimates with 
elevation angle are shown in Figure 1. The VAD technique has the highest 
efficiency of the techniques for all elevation angles. In addition, the VAD 
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maintains reasonable efficiency to larger elevation angles (75°) than either 
fixed beam or VED. 

Figure 1 shows efficiencies monotonically increasing as 6 gets smaller. 

But for measurements at a constant height the range increases when 0 
decreases. Because echo power reduces in proportion to the inverse square of 
range (assuming the echoing layer scatters isotropically and is horizontally 
homogeneous) the signal- to-noise ratio (SNR) falls as 6 decreases. If 
measurement errors are solely due to thermal noise and SNR is less than one, 
ZRNIC 1 (1979) shows that measurement variance o_ 2 is proportional to (SNR)” * 

The effect of decreasing SNR as 0 e becomes small is to increase o e ^ as esc 4 0 e ; 
the efficiency of the horizontal wind measurements thus vanishes as 0 e goes to 
zero. However, the variance o £ 2 includes meteorological effects such as 2 2 

turbulence that, in our experience, places a lower hound on of about 1ms 
Because we are mainly concerned in this paper with elevation~angles larger than 
40° and, consequently, ranges less than about 30 km, a. can be regarded as a 
constant. b 

The efficiencies of the vertical velocity estimates are shown in Figure 2. 
The VED has the highest efficiency, but for large elevation angles the VAD is 
comparable. The fixed beam has n a constant efficiency of since the number of 
vertical estimates is fixed at 

The biases of the estimates show a linear dependence on the height h. To 
normalize bias errors, we assume h = 1 km, so the bias for greater heights can 
be simply computed. The biases depend on the value of unknown spatial „ 1 
derivatives. Following WALDTEUFEL and CORBIN (1979) we use U x = v - l(T s a~ L 
and w x *= w = 10”^s”^ as maximum values. The biases thus computed y are shown in 
Figure 3 f£r the horizontal components and in Figure 4 for the vertical 
component. The asymmetry of the beam locations about the vertical for the fixed 
beam technique produces a horizontal wind bias due to u x and v y which decreases 
as cot0 e . The horizontal wind biases for the VAD and VED are the same and are 
constant with elevation angle. For the fixed beam technique the vertical velo- 
city is not biased by any derivatives. The vertical wind bias in the VED and 
VAD decreases with increasing elevation angle. 

In conclusion, the vertical velocity variance and the bias errors can be 
decreased by using larger elevation angles. The variance for horizontal 
components increases with elevation angle but can be controlled to an extent by 
data averaging. Because bias increases with height, the higher altitudes may 
require a vertical measurement for vertical velocity. 

(a) Example 

The bias the variance equations can be used to choose an elevation angle 
for profiling. For example, suppose we wish to prof ile the winds at 5 km using 
360 measurements with a £ » 1 m*s”l. The root mean square errors (bias squared 
plus variance) for the horizontal and vertical components are shown in Figures 5 
and 6 for each of the techniques. We have kept vertical and horizontal errors 
separate because vertical velocity is much smaller and requires greater ^ 

accuracy. I| we require horizontal and vertical velocity accuracies of 1 m-s~ 
and 0.1 m*s respectively, we would use an elevation angle between 83° and 85° 
for the VAD and 77° and 81° for the VED. Because of the bias error, the fixed 
beam horizontal wind error is 1.5 m- s” 1 or larger* 

POSSIBLE IMPROVEMENTS 

The analysis of the previous section suggests some simple improvements that 
can be made tcJ increase the accuracy of the measurements. 
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(a) Fixed Beam with Error Minimization 


The wind estimate efficiencies for the fixed beam technique can be improved 
slightly be collecting a specified number N of vertical data. Minimizing the 
first diagonal element of the matrix in (Al.l), which is VAR(u Q ) « VAR(v q ) 
gives 


, 7 sine 
N . m e 

n /2+sine 
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so N would vary from y 
efficiencies are shown 
that f or 0 e > 45°, the 
improved for w Q . 


( 12 ) 


at 0 e * 45° to about 0.4 n for e near 90°. The 
in Figure 7 for u ? and Figure 8 For w 0 . It can be seen 
efficiency of estimating u is unchanged but is slightly 



Figure 7. Horizontal wind estimator efficiency for fixed 
beam, fixed beam with horizontal error minimization, and 
fixed beam with equally spaced, off-vertical beams. 



Figure 8. Vertical wind estimator efficiencies for the 
same techniques as described in Figure 7 caption. 
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2 2 

The size of the variance contribution (i.e., cr £ tan N) from the 
vertical velocity bias removal appears to indicate that it might be better to 
ignore the bias error (w Q tan8 e ). However, for observation time intervals of 
several minutes, a meso scale value of w Q should be used* For large elevation 
angles, (0 e > 75°) the bias error could be several meters per second or larger. 
For tropospheric observation under all conditions, the bias should be removed if 
a horizontal velocity accuracy of 1 m*s 1 is required. 

(b) Three Off -Vertical Beams 

For some applications, ground clutter presents a problem for vertical 
measurements. The fixed beam technique with three off-vertical beams with 
elevation 8 e and azimuths 0°, 120°, 240° is analyzed in Appendix lb. The 
analysis shows that the variances (and efficiencies) of the estimates are 
identical to those for the VAD technique (compare Figures 1, 2 with Figures 7, 
8). The biases for the u 0 estimate is very similar to the fixed beam with one 
vertical but the w 0 estimate is biased as shown in Figure 9. 

(c) Application of the Continuity Equation to VAD Data 

Vertical winds as small as few centimeters per sec are important in fore- 
casting and, as noted earlier, w Q should be estimated with more accuracy than 
the horizontal components. Because of ground clutter it may become very 
difficult to estimate the radial component of air motion when the beam is 
pointed near the vertical since the radial velocities will have values close to 
zero. We now show that by assuming a linear wind field and applying the mass 
Continuity equation, we can estimate vertical wind, averaged over the circle of 
measurement, with the required accuracy. When mass continuity is applied, we 
will call the technique indirect whereas the previously discussed techniques 
(e.g. , VAD) are direct measurements of w. 

Two VAD modes in which vertical soundings can be made are fixed 8 e 
variable r, and variable 0 e fixed r. With variable r however, the horizontal 
area for which w Q is representative varies, so we prefer the second mode. 
Divergence is estimated by applying Gauss's theorem to the volume V (see Figure 
10) enclosed by the area at constant range from the radar and the area S 2 
at constant height (DOVIAK and ZRNIC^, 1983). Applying mass continuity and 
integrating gives an areal averaged w. 
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Figure 9. Vertical wind biases for the same techniques as 
described in Figure 7 caption. 
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Figure 10. Geometry to estimate vertical velocity 
averaged over the circular area S 2 . 
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where p « gH/RT is the average lapse rate of air density versus height and 
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is the average radial velocity around the circle of measurement. 

To fix the number of measurements at n, we assume M values of v r are 

made on each of L circles spaced at intervals Ah from h»0 to h=h • Then 
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If all the radial velocities are independent with the same uncertainty, then 
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For a direct measurement with a vertical beam VAR[w q ] « VAR[v r ]/N. If we 
require that, for our maximum height h^, VAR[w} « VAR[w q ] , then the range can 
be found by solving ° 
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Because of the accuracy needed for vertical velocity estimates, the number N of 
vertical data will be much larger than the n-N data for horizontal wind 
component estimation. Thus in (17) we can assume n/N~l. Using 1% ■ 10 km, 

T 0.113 km -1 in (17) gives r~40 km, so 0em~^°* For heights lower than h , 
VAR[w] is less than for a direct measurement. m 


To compare the variances for direct and indirect measurements, assume that 
we have n/10 measurements at each level for estimating u 0 , v Q , and w Q at each 
of 10 levels spaced 1 km apart. Solving (16) assuming a specified VAR[w] ■ 10“^ 
m 2 s~ 2 at 1% * 10 km with VAR[v r ] * 1 m s“ 2 gives n = 1080. So there are 108 
data at each level and, from (8) VAR[u 0 ]~2xl0~ 2 m 2 s~ 2 . With 108 data at each 
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level for the direct VAD measurement techniques, the VAR[w 0 ] is larger than 
that obtained using the indirect method even with 0 e = 90°. Because w Q needs 
to be estimated with better accuracy than Uq, v q , consider that most measure- 
ments are made with a vertical beam. We need at least 4 data spaced 90° apart 
on a circle at each level in order to remove the bias in Uq, v q due to v x 

(Equation 7). Thus from (8) we have VAR[a^]>0.5 m 2 s~ 2 which should be 
satisfactory for horizontal wind estimates but which is an order of magnitude 
larger than obtained by the indirect method. Furthermore, we have at most 104 
da ta^ available at each level for the vertical beam with the consequence that 
VAR[w 0 ]-10~ 2 nr • s” 2 , an order of magnitude or more (at h<10 km) larger than 
obtained using the indirect method. A comparison of these variances is shown in 
Figure 11. 

In reality, the variances for the indirect measurement technique may not be 
this much smaller, since we have neglected the dependence of signal strength 
(and VAR[v r ] with range. However, this analysis has shown that the indirect 
technique does not require a tradeoff between vertical and horizontal variance. 
It offers the advantages of low variances, an areal averaged vertical velocity, 
and requires no assumption about the spatial structure of the wind to measure 
vertical velocity. 

SUMMARY AND CONCLUSIONS 

We have examined the errors in three radar techniques (three fixed beams, 
VAD, and VED) used to directly measure the three components of the wind. 
Equations were derived for the bias and variance of the uniform wind components 
estimates under the. assumption of a spatially linear, time invariant wind field 



Figure 11. Comparison of wind estimate variances 
• for direct and indirect techniques. 
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and a constant radial velocity measurement error. The measurement errors 
produce variance in the estimates and the linear wind shear biases the 
estimates* The variance of the estimates can be reduced by averaging more 
measurements but the biases cannot. Thus, for these direct measurement 
techniques, the selection of an elevation angle for simultaneous observation 
requires a compromise based on the required accuracy of the measurement. 

We have also examined the errors for an indirect measurement technique 
based on Gauss* s theorem with an equation of continuity constraint* Two 
advantages this indirect technique offers are that it does not require any 
assumptions about the spatial structure of the wind to measure vertical velocity 
and that its error variance can be smaller because it does not require a 
vertical and horizontal measurement variance compromise* 
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APPENDIX 1» ANALYSIS OF THE FIXED BEAM TECHNIQUES 

For the fixed beam technique with one vertical, and two off-vertical beams at 
elevation 0 e , azimuths 0° and 90°, we find, following the notations of KOSCIELNY 
et al. (1982), that 

cos0 0 sin8 e 
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Since this is a symmetric matrix, we have not entered the identical terms below 
the diagonal elements. We invert P^Po by performing a sequence of row operations 
to reduce it to the identity matrix. Performing this same sequence on the identity 
matrix reduces it to [P 3 T P 3 ]~ 1 (ANTON, 1981). Thus 
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For an equal number of measurements on each beam, 3N=n and 
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The biasing of the estimates by the derivatives of the linear wind can be com- 
puted from the alias matrix of (3). Since the analysis is for constant height, all 
vertical derivatives cause no bias. Thus, for equal numbers of measurements. 
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T —1 T 

A 34 = [P nm P nm l P nin and performing the multiplication gives 
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For the fixed beam technique with three beams at elevation 0 and azimuths 
120°, and 240°, we find that 6 
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The alias matrix can be computed as before. For this case the predictor 
function for deformation is nonzero so 
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where a * /3/2 and b « 1/2. Computing A as before. 
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APPENDIX 2. ANALYSIS OF THE VAD TECHNIQUES. 

For the VAD technique, n radial velocity data are collected on a circle at 
height h. So 
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where all summations are for i~l, 2,..., n. To simplify P ? nia > we approximate 
the summations by integrals. For example, 
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A similar evaluation of the remaining summations gives 
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To find the bias caused by neglecting the parameters of a linear wind field, 
re compute the alias matrix where K 1 =[u + v,u,v,w,w]. For analysis 
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erf arming the matrix multiplications in (5), approximating summations by integrals. 
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APPENDIX 3. ANALYSIS OF VED TECHNIQUE. 

In the VED technique, the azimuth is fixed and the elevation angle is scanned. 
The range r is selected so the data are for a constant height h. To measure both 
horizontal components, another azimuth, preferably at 90° to the first, must be 
scanned. We assume y data are collected for each scan. For convenience, we 
introduce the horizontal distance s which is directed along the azimuth 4> • We 
make estimates of a horizontal component u G (directed along s) and the vertical 
component w q . The predictor function matrix is 



where the summations are for 1=1, 2,,.., — and 9^ is the elevation angle. 
Approximating the summations by integrals, 
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The variance of w is halved since we assume the results from the two scans 
will be averaged. 

The bias by the linear terms is again computed by the alias matrix. The 
prediction functions corresponding to the excluded parameters are 
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Performing the matrix multiplications, approximating summation by integrals and 
using h~r sin0 gives 
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